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We extract the magnetic polarizability from the quadratic response of a hadron's mass shift in progressively 
small static magnetic fields. The calculation is done on a 24 x 12 x 12 x 24 lattice at a = 0.17 fm with an 
improved gauge action and the clover quark action. The results are compared to those from experiments and 
models where available. 
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Polarizabilities are important fundamental 
properties of particles. They determine dynami- 
cal response of a bound system to external per- 
turbations, and provide valuable insight into in- 
ternal strong interaction structure. We discuss a 
lattice calculation for static magnetic polarizabil- 
ity P in the quenched approximation. This work 
is an extension of a lattice calculation for the elec- 
tric polarizability [0. For an updated calculation 
of the electric polarizability, see Q . 

For small external magnetic fields, P can be 
defined via the energy (mass) shift of a particle 
in the fields 
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We determine the magnetic polarizability by cal- 
culating the mass shift of hadrons 



Am = m{B) — m(0) 



(2) 



where m(B) is the hadron mass in the presence 
of an magnetic field B. We fit the data with the 
polynomial 

Am{B) = ciB + icaS^ (3) 

*presented by L. Zhou at Lattice 2002. 
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We calculate mass shifts both in the field B and 
its inverse —B, then average them. The magnetic 
polarizabilities are the negative quadratic coeffi- 
cients 



(4) 



The gauge action is the zero-loop, tadpole- 
improved Liischer-Weiz action with a=0.17 fm (or 
l/a=1159 MeV) set from the string tension. The 
tadpole factor is determined by the average pla- 
quette as uq ~ 0.877. The clover quark action was 
used to put fermions on the lattice. 
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where Sw is the tadpole-improved Wilson 
fermion action. 

In order to place an magnetic field on the lat- 
tice, we construct an analogy to the continuum 
case. There, the fermion action is modified by the 
minimal coupling prescription ^ + iq^n, 
where q is the charge of the fermion field and 
is the vector potential describing the external 
field. On the lattice, the prescription amounts 
to a phase factor e*'"''^f to the link variables [Q. 
Choosing A2 = Bxi and A^) = Ai = — Q 
a constant magnetic field F12 = diA2 — d2Ai = 
B3 = B can be introduced in the z-dircction. We 
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Figure 1. Mass shifts of the neutron as a function 
of the magnetic field B at the six quark masses. 

introduce two dimensionless parameters: one is 



7/ = qBa 



(6) 



where q = Qe with Q = ±1/3, ±2/3. The other 
is the integer lattice length, 



p — xi/a. 



(7) 



In terms of these two parameters the phase factor 
becomes: 



JaqA2 



(1 + i?7p). 



(8) 



where we have linearized the phase factor to 
mimic the continuum prescription. Therefore, 
the field strength we use will be chosen to satisfy 
the linearization requirement. To summarize, our 
method to place the external magnetic field on the 
lattice is to multiply each gauge field link variable 
U^{x) in the z-direction with a a;-dependent fac- 
tor: 



U^{x)^(l + i7^p)U^{x) 



(9) 



This factor is applied only to the Sw part of the 
action, the clover term is untouched. 

Periodic boundary conditions were imposed on 
the fermion fields in y and z directions. In the 
X and t directions, we have employed fixed (or 
Dirichlet) boundary conditions, i.e. the fermion 
fields are not allowed to propagate across the 
edges of x and t. To minimize the non-periodic 
boundary effects, we use a lattice size of 24 x 
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Figure 2. Magnetic polarizability of the neutron 
and proton as a function of in physical units. 

12 X 12 X 24 with the fermion source location 
{x,y,z,t) = (12,1,1,3). We analyzed 100 con- 
figurations. 

Fermion propagators were constructed at 
six K values, k=0.1219, 0.1214, 0.1209, 0.1201, 
0.1194, 0.1182 corresponding to m^=450 to 820 
MeV. The critical value Kc = 0.1232(1). We 
use six different values of the parameter 77 = 
0.0, +0.00036, -0.00072, +0.00144, -0.00288, 
+0.00576. The 77 values in this sequence are re- 
lated by a factor of —2. Thus we are able to 
study the response of a hadron composed of both 
up and down (strange) quarks to four different 
nonzero magnetic fields. 

Typical mass shifts are displayed in Fig. || in 
the case of neutron. There is good parabolic 
behavior going through the origin, an indication 
that contamination from the linear term has been 
effectively eliminated by averaging results from B 
and —B. The factor 



0.358 X 10 



«cm3 



(10) 



was used to convert the lattice numbers into phys- 
ical polarizability in units of lO^^fm"^. The curves 
in Fig. ^ and Fig. || represent our preliminary re- 
sults for neutron, proton, tt^ , and K^. They are 
smooth functions of the quark mass. We did not 
attempt an extrapolation in the quark mass. 
In Table we show some experimental and 
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Figure 3. Magnetic polarizability for pion and 
kaon. 



theoretical values for the magnetic polarizabili- 
ties f3 of neutron, proton, pion and kaon. The 
/3 for neutron is not well-determined. A recent 
experimental value ||^ is 2.7 with a big error bar, 
while an earlier experimental data analysis has 
it at 11. The range from models and experiments 
is from 1.22 to 14. Our lattice simulation result 
is at the higher end of this range. For the proton, 
the range is 1.6 to 4.4, while our lattice result lies 
in the range. For tt", there is scant data. The 
number from Chiral Perturbation Theory calcu- 
lation to 0{p^) is 0.5 to 1.7. Our result is quite 
high, about 20. For K'^, there is one theoretical 
value —6.2 from the NJL model. Our result is 
around -20. 

In summary, we have computed the magnetic 
polarizability of hadrons for the first time on the 
lattice using the external field method. In addi- 
ton to the particles described here, we have com- 
puted (3 for other mesons, the other members of 
the octet baryon, and the entire decuplet baryon 
family. We plan to present the final results else- 
where. Also under way is the extraction of mag- 
netic moments from the same data set. 
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